function P3_Gene_Regulatory_Network
close all;clear all;
% reproduce the figure 4.8c
% show three equilibrium points
GeneRegulation_Nullcline_with_Jacobian(3,2);

% superimpose the direction fields on the same axes
Direct_Fields(3,2);

% plot three trajectories
% starting from (0,0)
Trajectories(3,2,0,0);

% starting from (0.5,0)
Trajectories(3,2,0.5,0);

% starting from (0,0.5)
Trajectories(3,2,0,0.5);

end

